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Abstract 

We perform a study of the Standard Model (SM) fit to the mixing quantities AMb^, and 
ATbs/AMb^ in order to bound contributions of New Physics to Bs mixing. We then use this 
to explore the branching fraction of — )• in certain models of New Physics (NP). In most 

cases, this constrains NP amplitudes for Bg — ?■ fJ-^fJ^^ to lie below the SM component. 
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I. INTRODUCTION 



We report here on a study of New Physics (NP) predictions for Bg — ?■ n'^^~. The 
Standard Model (SM) prediction for Bg — ?■ fj^^fi~ is currently smaller than the experimental 
branching fraction limit {l| of by about a factor of 15. This presents a window of 

opportunity for observing New Physics (NP) effects in this mode. 

This topic is particularly timely in view of experimental indications of NP effects in both 
the exclusive decay 5^ — )■ + $ [2] (for recent CDF results, also see Ref. pj]) as well as 



the inclusive like-sign dimuon asymmetry observed in pp — )■ jjjj, + X ^|. Moreover, future 
work at LHC-B, e+e^ Super B-factories and ongoing CDF & DO measurements at Fermilab 
(see the discussion following Eq. (E])) is expected to markedly improve the current branching 
fraction bound. 

n 

Our strategy in this paper is somewhat reminiscent of our recent study [5] noting that 
in some NP models the mixing and — ?■ fj^~^fi~ decay amplitudes have a common 
dependence on the NP parameters. If so, one can predict the yU+yU" branching 

fraction in terms of the observed AM^ provided that much or all of the mixing is attributed 
to NP. This is a viable possibility for mixing because the Standard Model (SM) signal 
has large theoretical uncertainties and because many NP models can produce the observed 



mixing jo]. 

For AMb, the situation is very different. Here, the SM prediction is in accord with the 
observed value {e.g. see Refs. [7|, |8| and papers cited therein). In fact, the analysis described 
below (c/. see Eqs. (inD,(II3])) gives |AM|j^^Va^b!^^I < 0.20, which demonstrates just 
how well the SM prediction agrees with the experimental value of AMb^. In view of this, 
our SM expression for AMb^ will be given at NLO [ol, [lo| whereas LO results will suffice 
for NP models. As regards the corresponding width difference AF^^, the experimental and 
theoretical uncertainties are still rather significant {viz Sect. II-C). 

In those NP models where mixing and Bs — ?■ fJ^^fJ^' arise from a common set of parameters, 
the severe constraint on any NP signal to Bs mixing places strong bounds on its contribution 
to Bb^^^+^--^ In fact, we shall find the constraint can be so strong that for some NP models 
the predicted Bg — )■ /x^/x^ branching fraction lies well below the SM prediction. 



^ In particular, Ref. [7j considers the possibility, not covered here, on effects of so-called minimal flavor 
violation which affect the quark mixing-matrix elements. 
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The first step in our study (c/ Section II) will be to revisit the SM predictions for mixing 
in the 6-quark system by using up-to-date inputs. We carry this out for the two mixing 
quantities AM^^ and AFg^/AM^^. The former in turn yields phenomenological bounds on 
NP mixing contributions which in certain models can be used to bound the magnitude of 
the Bs —7- fi^fi^ decay mode. We also update the SM branching fraction for Bg by 
using the observed Bs mixing as input. Then, in Section III we discuss general properties of 
NP models with tree-level amplitudes. In Section IV, we explore various NP models such as 
extra Z' bosons, family symmetry, R-parity violating supersymmetry, flavor-changing Higgs 
models, and models with the fourth sequential generation. Our concluding remarks appear 
in Section V, and some technical details are relegated to the Appendix. 

II. UPDATE OF Bs MIXING AND Bs IN THE STANDARD MODEL 

We begin by considering the SM predictions for Bg mixing. This step is crucial to obtaining 
bounds on NP contributions. We also use the Bs mixing signal as input to a determination 
of the branching fraction for Bs — ?■ fi^fi". 

A. Inputs to the Analysis 

The work in this Section takes advantage of recent progress made in determining several 
quantities used in the analysis. We summarize our numerical inputs in Table I, along with 
corresponding references. Included in Table I is an updated determination of the top quark 



Mb, = 5366.3 ± 0.6 MeV [Ij 


TB, = (1.425 ± 0.041) X 10-12 g py 


AMb, = (117.0 ± 0.8) X 10-13 


ATBjTB,=0m2tlf,lm 


XB^ = 0.776 ±0.008 yy 


XB, = 26.2 ±0.5 [1] 


^(pole) ^ ^^3_^ ^ ^ 3 ^^^^ 


as{Mz) = 0.1184 ± 0.0007 [Jj 


Jb, = 0.2388 ± 0.0095 GeV \13\ 


Jbs^Bb, = 275 ±13 MeV [13] 


|V,.|= 0.0403l°X [1] 


m = 0.999152^0;— [1] 



TABLE I: List of Input Parameters 
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pole mass 



11 1 m[^"^^^ which in turn is used to determine the corresponding running mass 
TTT'ti'fnt) IJl along with several decay constants and B-factors as evaluated in lattice QCD. 
For definiteness, we have used values appearing in Ref. jisl. This area is, however, constantly 
evolving and one anticipates further developments in the near future jl5i]. Our values for the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |Vts| and \ Vih \ a re taken from Ref. 
Similar values occur for the global fits cited elsewhere [e.g. Refs. [IgI. \\\)- 



B. /\Mbs 
The PDG value for AM^^, 

^^Jexpt) ^ ^^^^ Q ^ Q ^ ^Q_13 ^^y ^ j-^^j 

is a very accurate one - the uncertainty amounts to about 0.7%. The NLO SM formula, 

is arrived at from an operator product expansion of the mixing hamiltonian. The short- 
distance dependence in the Wilson coefficient appears in the ^ scale-insensitive combination 
VBsSoi^t), where the factor So{xt) is an Inami-Lin function [18] (with Xt = m^{rht)/M^) 
and ffit{fht) is the running top-quark mass parameter in MS renormalization. In particular, 
we have fht{rht) = (163.4 ± 1.2) GeV which leads to So{xt) = 2.319 ± 0.028. Using the same 
matching scale, we obtain rjBs = 0.5525 ± 0.0007 for the NLO QCD factor. 
Our evaluation for AM^^'^ then gives 

AMi^^ = (125.2111?) X 10-13 , (3) 

which is in accord with the experimental value of Eq. ([1]). The theoretical uncertainty in 
the SM prediction of Eq. ([3]) is roughly a factor of sixteen larger than the experimental 
uncertainty of Eq. ([2]). The largest source of error occurs in the nonperturbative factor 
BbsIbs^ followed by that in the CKM matrix element Vts- The asymmetry in the upper and 
lower uncertainties in AM^j^'^^ arises from the corresponding asymmetry in the value of Vts 
cited in Ref. 

Finally, we note in passing that for the ratio AMb^/ AMb^ the experimental value is 
0.02852 ± 0.00034 whereas the SM determination gives 0.02835 ± 0.00187. This good agree- 
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ment is not surprising since the ratio AMb^/AMb^ contains less theoretical uncertainty 
than AMbj^ or AMb^ separately. 



C. The Ratio ATbJ^Mb, 

The above discussion of AM^^'' sets the stage for analyzing NP contributions to Bg — )■ 
/i"*"/!". There is, in principle, a second approach which instead utilizes AF^^. The PDG 
value for the 5, width difference is AT^^T^ = OM2tom7 x lO^^g-i^ Together with Eq. (P, 
this gives^ 

^(expt) ^ Bs ^ -omi ^ 34 g ^ 20.0 X 10"^ . 4 

^^^(cxpt) (17.77 ± 0.12) X 1012 s-1 ^ ' 

whereas the corresponding SM prediction from Ref. ^] is r^^^'' = (49.7 ± 9.4) x 10~^. In 
contrast to the mass splitting AMb^, the theoretical uncertainty in the ratio ATbJAMb^ 
is much smaller than in the current experimental determination. Nonetheless, this situation 
is expected to change once LHCb gathers sufficient data. As such, we would expect a highly 
accurate value of APg^''*^ to eventually become available. We propose that it could be 
applied to the kind of analysis used in this paper as follows. We define a kind of mass 
difference VMb^ as 



AM 



(thy) 



^^B. - Arg-*) . (5) 



Ar 



(thy) 
Bs 

(thy) 



The point is that if NP contributions are neglected in AB = 1 transitions, then AF^ 
is purely a SM effect. In addition, the ratio AM^^^/AT^^^^ will be less dependent on 
hadronic parameters than either factor separately. 

This quantity is also important in the scenarios where NP contributes a significant CP- 
violating phase to AMb^- In this situation, AFg^^*"* will be reduced compared to its SM 
value AF^^^'* by a factor of cos 2£, where ^ is related to the relative phase between the SM 



and NP contributions to AMb 



19|. 



At the very least, the relation in Eq. ([5]) would be of interest to analyze the NP issue 
using both quantities AMb^ and the above VMb^- 



Using instead the recent CDF evaluation AF^^ ^ = 0.075 ± 0.035 ± 0.01 x lO^^ g-i implies A'^'^vt) 
(42.2 ± 20.5) X 10-", consistent with the value in Eq. (g]). 
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D. B., 



PDG entries for Bb^ 



are 



(expt) 



< 4.7 X 10- 



and 



B 



(expt) 



< 2.8 X 10' 



(6) 



with no experimental limit currently for the Bg — t- r+r transition. Data collected by the 
DO and CDF collaborations will improve the above brancing fraction limit. For example, 

?(D0) 



the DO collaboration reports < 5.1 x 10' 



with an anticipated limit of eleven 
□ 

times the SM prediction and similarly for the CDF collaboration [20|. 



Since the long distance (LD) estimate for the branching fraction of Bg — ?■ fi^fi in the 



SM gives B 



(LD) 

Bs^^l+^l' 



6 X 10^1^ 



2l[ |. we consider only the short distance (SD) component 



in the following. Using Eq. ([2]) as input to the SD-dominated Bg — )■ transition (see 

also Ref. {?]) we arrive at 



B 



(SM) 

Bs-^fi+fi- 



^Mb, tb 



Ar]B,BBsT^^ 



ml 

1-4- 



'Ml 



1/2 



Soixt) 



(7) 



where Y{xt) is another Inami-Lin function 18|]. Expressing B 



+ - in this manner serves 



to remove some of the inherent model dependence. Numerical evaluation gives 



B 



(SM) 



(3.33 ±0.21) X 10" 



(8) 



The major sources of uncertainty, ordered by magnitude, arise from the factors Bb,, then 
the lifetime tb^ and finally the top-quark mass value. 

III. STUDY OF NEW PHYSICS MODELS 

In this section, we first obtain a numerical (Icr) bound on any possible New Physics 
contribution to AMb^ ■ We then use this to constrain couplings in a variety of NP models 
and thereby learn something about the Bg — t- fi^fi'^ transition. 



A. Constraints on NP Models from Bg Mixing 



As shown in Ref. 22], New Physics in AB = 1 interactions can in principle markedly 
affect AP^. The logic is similar to that used in Ref. [231 regarding the possible impact of 
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NP on ATd. Since, however, in Bg mixing such models are not easy to come up with, one 
can simply assume that AB = 1 processes are dominated by the SM interactions. Thus we 
can write 



AMs^ = AMg^) + AmS'P cos <f> , (9) 

If the AB = 1 sector were to contain significant NP contributions, then the above relation 
would no longer be valid due to interference between the SM and NP components. 

As can be seen from Eq. interference between the SM and NP components may also 
occur in the presence of a CP- violating phase in the NP part of the mixing amplitude 24 1. 



This large NP phase could markedly affect AP^^^*'' even in the absence of a NP contribution 
to the on-shell AB = 1 transitions (recall that AP^^''*^ depends explicitly on the cosine of 



19|) 



the CP- violating phase ^ [sl, [3]; the explicit relation between and ^ can be found in 
It is therefore more reasonable to use AP^^^*'' in studying those scenarios with a large NP 
phase. The appropriate strategy here would be to use AP^^^*'' and AP^^'^^ to extract the 
phase ^, eliminate cos0 from Eq. iQ, and then extract AM]j^^^ in order to relate it to 
the rare leptonic decay rate. To do so, however, will require a significant reduction in the 
experimental uncertainty of AP^^^*'' . Alternatively, CP- violating phases could be extracted 



at LHCb from the studies of Bg — )■ J /ipip transition 2^. We shall defer those studies to a 



future publication 25|. In this paper we shall assume that the phase in the NP component 



of AMb^ is sufficiently small (although not necessarily negligible), 

AMb^ = AMg^) + AM^^f^ . (10) 
Accounting for NP as an additive contribution, 

^^(expt) _ ^^(SM) ^ ^^(NP) ^ ^^^^ 

we have from Eqs. ([I]),(l3]), 

AM^f ) = (-8.2^^3-8) X 10-13 . (12) 

The error in AMi^^P*) has been included, but it is so small compared to the theoretical error 
in AM^^^^ as to be negligible. The la range for the NP contribution is thus 

AM^^f^ = (-20.9 +5.6) X lO^^^ GeV . (13) 



To proceed further without ambiguity, we would need to know the relative phase between 
the SM and NP components. Lacking this, we employ the absolute value of the largest 
possible number, 

lAM^sf^l < 20.9 X 10-^^ GeV , (14) 
to constrain the NP parameters. 



B. Generic NP Models with tree-level amplitudes 

New Physics can affect both Bg mixing and rare decays like Bg — t- fj^^fi^ by engaging in 
these two transitions at tree level. In this section we will, for generality, consider a generic 
spin-1 boson V or a spin-0 boson S with flavor-changing and flavor- conserving neutral current 
interactions that couple both to quarks and leptons. The bosons V and S can be of either 
parity. This situation is frequently realized, as in the interactions of a heavy Z' boson or in 
multi-Higgs doublet models without natural flavor conservation. 

Spin-1 Boson V: Assuming that the spin-1 particle V has flavor-changing couplings, the 
most general Lagrangian can be written as^ 

•Hy = g'vAl/LV' + g'yAl/nV^ + gvihl.SLV'' + gv2bRi,SRV^ + h.c. . (15) 

Here is the vector field and the flavor of the lepton i' might or might not coincide with 
i. It is not important whether the field V"^ corresponds to an abelian or non-abelian gauge 
symmetry group. Using methods similar to those in Ref. Q], we obtain 



^= 3M? 



2 /\/r_ r i; J 

Ci{fi)B, + Ce{fi)Be - -C2(/i)52 + -C3{fi)Bs 



V 

(V) 



(16) 



where the superscript on AM^J denotes propagation of a vector boson in the tree amplitude. 
The Wilson coefficients evaluated at a scale fi are related to the couphngs gyi and gv2 as 



Ci(/i) = r(/i. My) gl, , CM = I [K/U, Myf'^ - r(/i, My)-\ 

C^ifi) = 2 r{fi, Myy/^gy,gy2 , CM=r{fi,My) g^y^ , 



gvigv2 



^ Throughout, our convention for defining chiral projections for a field q{x) will be qL,R{x) = (l±75)(?(x)/2. 



8 



where (presuming that M > rrit and /i > nib), 



r{fi,M) 



6/23 



(17) 



Similar calculations can be performed for the — £^£^ decay. The effective Hamiltonian 



in this case is 



H 



(V) ^ _J_ 



(18) 



gvig'viQi + 9vi9v2Q7 + g'vi9v2Q2 + gv2g'v2Q&\ > 

where the operators {Qi] can be read off from those in Ref. with the label changes c — )■ s 
and u ^ b. This leads to the branching fraction, 



(V) 



Afflf . 1 2 1 / ' 1 2 

-L - \gvi - gv2\ \gvi - gv2\ 



(19) 



Clearly, Eqs. ( !T6|) .( !T9|) can be related to each other only for a specific set of NP models. 

Spin-0 Boson S: Analogous procedures can be followed if now the FCNC is generated 
by quarks interacting with spin-0 particles. Again, the most general Hamiltonian can be 
written as 



= g'si^dnS + g's2^dLS + gsihsnS + gs2hRSLS + h.c. . 
Evaluation of AMg ^ at scale // = gives 



(20) 



AM, 



(s) 



24M| 



Tie 



(21) 



with the Wilson coefficients defined as 



CM 

C4(/i) 

CM 
CM 

CM 



-2r{fi,Ms) ^ gsigs2 = C-M gsigs2 



r_(/i,Ms 



gk = CM gs2 



1 



8^241 

'1 



[r+(/i, Ms) - r^fi, Ms)] gl, = CM gk 



(22) 



^= r+(u, Ms) + - + ^= r_(u, Ms) 



%i = Crin) gsi 



8^241 



[r+(/i, Ms) - r_(/x, M^)] (7^1 = Cs{fi) g^ 
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where for notational simplicity we have defined r± = 7-(i±v^24T)/6_ jsjQ^g thai Eq. f l2T|) is true 
only for the real spin-0 field S. If S* is a complex field, then only operator Qs will contribute 
to Eq. (ED). 

The effective Hamiltonian for the 5° — decay via a heavy scalar S with FCNC 
interactions is then 



9sig'siQ9 + gsig's2Qs + g'si9s2Q3 + 9s2g's2Q^ > (23) 



and from this, it follows that the branching fraction is 

j(S) _ Pb^b, 



B 



BO- 



1 - 



Ami 



— \9si — gs2\ 



2 ' 



^-jjT-j + \g'si-g's2\' 



(24) 



Note that if the spin-0 particle 5* only has scalar FCNC couplings, i.e. gsi = gs2, no contri- 
bution to 5° — £^£^ branching ratio is generated at tree level; the non-zero contribution to 
rare decays is instead produced at one-loop level. This follows from the pseudoscalar natuTe 
of the i?<j-meson. 

Let us now consider specific models where the correlations between the Bs — Bg mixing 
rates and (in particular) the B^ — )■ rare decay can be found. 

C. Z' Boson 

Bs Mixing: The Bg mixing arising from the Z' pole diagram has the same form as in 



mixm 



AMb^ - - ^ , (25) 

where ri{mb, Mz') is a QCD factor which we take to be 

ri{mb,Mz') ~ 0.79 . (26) 

This is a compromise between ri(m5, 1 TeV) = 0.798 and ri(mfe,2 TeV) = 0.783. Solving 
for the Z' parameters, we have 

^ Z'sb I tSs I ^ Q no XX 1 n^ii 



Ml MBJlBBM^b,Mz' 



< 3.02 X 10-^1 GeV~2 (27) 
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upon using the constraint from Bs mixing. 

Bs — 7- Decay: This has aheady been calculated for — t- yU^/i^ decay in Ref. jsl- 

Inserting obvious modifications for Bg, we have from the branching fraction relation 

Eq. (39) of Ref. fl, 



B 



(Z') 



1 



B.^,+,' 16x/27rrB, ^ ^ Ml Ml Ml ' 
Upon inserting numbers, we obtain 



(28) 



B 



(Z') 



< 0.30 X 10^ 



1 TeVy 
Mz' ) ■ 



(29) 



This value is already below the corresponding SM prediction {B^B^fi+fi- ~ ^-^S x 10 ^) even 
if we take a Z' mass as light as Mz' — 1 TeV. 



D. R Parity Violating Supersymmetry 

One of the models of New Physics that has a rich flavor phenomenology is R-parity 
violating (RPV) SUSY. The crucial difference between studies of RPV SUSY contributions 
to phenomenology of the up-quark (see Q|) and down- type quark sectors is the possibility of 



tree- level diagrams contributing to i?s-mixing^ and Bs i^i decays 26|-|29| . If one allows 



for R-parity violation, the following terms should be added to the superpotential. 



Wf, = \\,uU^El + KjuUQ.Dl + \>:i^,UtDpi. (30) 



Here Q and L denote SU{2) l doublet quark and lepton superfields, and U , D and E stand for 
the SU{2)l singlet up-quark, down-quark and charged lepton superfields. Also, fc} = 
1,2,3 are generation indices. We shall require baryon number symmetry by setting A" to 
zero. Also, we will assume CP-conservation, so all couphngs Xijk and X'^^^ are treated as 
real. 

— 5^ Mixing: Neglecting the baryon-number violating contribution, the Lagrangian 
describing RPV SUSY contribution to B^ — B^ mixing can be written as 

= -K23^iLbRSL - K32^iLSBbL + h.C. , (31) 



^ We assume that there is no strong hierarchy between the RPV SUSY couplings that favors possible box 
diagrams. 
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where i = 1, 2, 3 is a generational index for the sneutrino. Matching to Eq. ( 120|) imphes 
that the only non-zero contribution comes from the operator Qs- Taking into account 
renormalization group running, we obtain for AM,, from the R-parity violating terms, 



AM 



m 



where M^^ denotes the mass of the sneutrino of zth generation and the function 



i^(C3,53 



^C3(/X,M,Ji?3, 



(32) 



(33) 



and the B-factor is defined in 



is defined in terms of reduced Wilson coefficient of Eq. 
Table [Til of the Appendix. 

Bg —7- /i^/U^ Decay: In RPV-SUSY, the underlying transition for Bg — )■ fi^fi^ is s + 6 — )■ 
/i+ + fi~ via tree-level w-squark or sneutrino exchange. In order to relate the rare decay 
to the mass difference contribution from RPV SUSY AM^ , we need to assume that the 
up-squark contribution is negligible. This can be achieved in models where sneutrinos are 
much lighter than the up-type squarks, which are phenomenologically viable. Employing 
this assumption leads to the predicted branching fraction 

2 / 



64 TT Tb, 



X 



E 



^i22^j32 



2 



Ml 



\ 



r_ 

Ml 



Ml 



+ 



E 



-^j22A-23 



Ml 



(34) 



In order to relate Bg — > yU+yU to AM^ in the framework of RPV SUSY, we need to make 
additional assumptions. In particular, we shall assume that the sum is dominated by a single 



sneutrino state, which we shall denote by z/^. In addition, we will assume that 



fc23 



A:32) 



which will reduce the number of unknown parameters. This assumption is not needed, 
however, if one wishes to set a bound on a combination of coupling constants directly from 
the experimental bound on -6^^^^+^-. Then, neglecting CP-violation, 

2 / 



f Bs '^Bs 



•^i22A-32\ / Mb, 



Ml, 



\ 



Ami 
Ml' 



(35) 



647r Tb, \ M? ; V m, 
where A; = 2 if an assumption that A'^23 = -^132 is made, and k = 1 otherwise. 

Since no Bg — t- signal has yet been seen, we can use the experimental bound to 

obtain an updated constraint on the RPV couplings. 



Afc22Afc32 < 5.5 X 10" 



M, 



100 GeV 



(36) 
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S.xlO- 



6.x 10"* 



a. 

T 



4.x 10"'* - 



2.X10- 




0.00 



0.01 



0.02 



0.03 



0.04 



FIG. 1: Branching ratio of B^o^^+fj,- as a function of RPV leptonic coupling Afe22 and sneutrino 
mass Mi^. = 100 GeV, 150 GeV, and 200 GeV (solid, dashed, and dash-dotted lines). The yellow 
shaded area represents excluded parameter space. 



Now, assuming = -^^32) can relate the branching ratio Bb, 

2 / 



B 



(W 



Ml 



M, 



Bs 



1 



Ml, 



\ 



j.m 

Ml Ml ■ 



(37) 



207r F(C3, V m. 

It is possible to plot the dependence of Bs^^fM+fi- on Afc22 for different values of M^;., which 
we present in Fig. [1] 



E. Family (Horizontal) Symmetries 



The gauge sector in the Standard Model has a large global symmetry which is broken by 



the Higgs interaction 



30|. By enlarging the Higgs sector, some subgroup of this symmetry 



can be imposed on the full SM lagrangian and the symmetry can be broken spontaneously. 



This family symmetry can be global 



3l| as well as gauged [32| . If the new gauge couplings 



are very weak or the gauge boson masses are large, the difference between a gauged or 
global symmetry is rather difficult to distinguish in practice 33|]. In general there would be 
FCNC effects from both the gauge and scalar sectors. Here we study the gauge contribution. 
Consider the family gauge symmetry group SU{3)g acting on the three left-handed families. 
Spontaneous symmetry breaking renders all the gauge bosons massive. If the SU(3) is 
broken first to SU(2) before being completely broken, we may have an effective 'low' energy 
symmetry SU{2)g- This means that the gauge bosons G = {Gi} {i = 1, ... ,3) are much 
lighter than the {Gk} {k = 4, . . . ,8). For simplicity we assume that after symmetry breaking 
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the gauge boson mass matrix is diagonal to a good approximation. If so, the hght gauge 
bosons G are mass eigenstates with neghgible mixing. 
The LH doublets 

(A 



(38) 



L ^ ' L ^ ' L 

transform as /g = 1/2 under SU{2)g, as do the lepton doublets 



(A 



(39) 



L ^ ' L ^ ' L 

and the right-handed fermions are singlets under SU{2)g- In the above, the superscript 'o' 

refers to the fact that these are weak eigenstates and not mass eigenstates. The couplings 
of fermions to the light family gauge bosons G is given by 



where / denotes the coupling strength and r are the generators of SU{2)g 
The fermion mass eigenstates are given by, first for quarks. 



(40) 



'A 




(A 


s 


= 






L 







( \ 




( A 




u 




and 


c 










L 


V J 



(41) 



and then for leptons. 



/ \ 




( A 




( \ 




( A 


e 












= u, 




and 








v) 


L 


V ) 


L 


vA 


L 





(42) 



The four matrices C/^, U(, and are unknown, except for 

UlUa = VcKM and UlUi = Vmnsp 



(43) 



where Vmnsp is the Maki-Nakagawa-Sakata-Pontcorvo lepton mixing matrix. The couplings 
of the gauge bosons relevant for the Bg system in the mass basis are: 



/ 
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[UblU*^bL'J^,SL + UslUl^SLlpbh - Ub2U*2bL'JfMSL - Us2Ub2SL bL 



The contribution to — B^ mixing is given by 



AM; 



(FS) _ '^MBjlBBr{mB, 



A 



C 



B 



2 ~^ 2 ~^ 2 
mi m^ m2 



(44) 



(45) 



where 



A = Re 
B = -Re 



(Ub, u:, + Ub2 u:,] 



{Ubi u:, - Ub2 u:,] 



(46) 



C = Re 



iUbiu:^ - Ub2u, 



s2J 



In a simple scheme of symmetry breaking 3J], one obtains mi = m^ and the square 
bracket in Eq. fH5]) becomes 



2 



5 

ml 



Although the matrices Ui {i = d, u, €) in principle are unknown, it has been argued that a 



reasonable ansatz 



35 1, which is incorporated in many models is Uu = I, U\ = Vckm- In 



this case^ one can simplify A, B and C further: 



< C ~ 1.6 X 10" 



Thus the Bg mixing becomes 



, fFS) 2MbJI Bb r(mb,M) P , 
AMlf^^ ~ ^■"^^ ^ ^ ' — - At 1.6 X 10"^ 

Bs 3 



2 



(48) 



(49) 



so that, substituting experimental bound AM^^p = AM'^ 



r(NP) 



2 

mf 



< 



(50) 



2MbsPb BBj{mb, M)1.6 X 10-^ 

The same above ansatz also implies that Uj = Umnsp and = I- Then the couphng of the 
gauge bosons to muon pairs is given by 



f 



+ * (-t^Mi u;, + u;,u,2) + (f/^i u;, - u,2 u;,) a, 



(J'L'-fXf^L 



(51) 



^ Here, we use values listed in Ref. [if. 
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The branching ratio for Bs — t- /i^/i is given by 



t^Bs^t^+fJ. 00_Tn / 



327rr 



ml 



mn 



Next we employ the approximation 
Utbm is the tri-bi-maximal matrix 



well-supported empirically) that f/^ 
'36! ] ■ Then Eq. ( 15T|) becomes 



(52) 



f^TBM, where 



-/ 



3^6^ 



(53) 



With this, the contribution to the branching ratio for Bg — )■ becomes 



Bs^tJ.+ li~ 



327rrB. 
MsJlmlf' 1.4x10- 



(1.1 X 10" 



X 0.04 



327rr7 



(54) 



The dependence on unknown factors in Eq. ( l54l) (z.e. [f /miY) can be entirely removed by 
using the bound in Eq. fl50|) to yield 



i3 



(FS) 

Bs^fj.+H 



< 



3.85ml 



ttMbsTbs {fBsBBsr{mb,mi) 
From the bounds of Eqs. 0121) . 0131) . we obtain 



2 \^Mb, 



(NP)|2 



(55) 



^R%^+„- < 0-92 X 10-^2 . 



(56) 



F. FCNC Higgs interactions 



Many extensions of the Standard Model contain multiple scalar doublets, which increases 
the possibility of FCNC mediated by flavor non-diagonal interactions of neutral components. 
While many ideas exist on how to suppress those interactions (see, e.g. j38-40|), the ultimate 
test of those ideas would involve direct observation of scalar-mediated FCNC. 

Consider a generic Yukawa interaction consisting of a set of Higgs doublets H„ {n = 
2, .., A^) with SM fermions. 



T^Y — ^UnQLiURjHn + X^^j^Q Li^ Bj Hn + Xfj^LliERjlln + h.C. 



(57) 
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where Hn = io'2H* and Qii (Lu) are respectively the left-handed weak doublets of an ith- 
generation of quarks (leptons). Restricting the discussion to Bg Mixing and — t- 
decay, we find that Eq. (157|) reduces to 

K = >^23nSLbR^l + Ai„6LS,J$° + Xi2npLf^R^n + h-C, (58) 

where $° = (0° + ) / -\/2. Bringing this to the form of Eq. f l2U]) and confining the discus- 
sion only to the contribution of the lightest 0° and a° states, we obtain 

- hRSLO^ + bLSRa'^ + i-^ 7^L/i/?a° + - + h.c. , (59) 

where ellipses stand for the terms containing heavier 0° and a° states whose contributions 
to AMb^ and Bb^^^+^- will be suppressed. 

If the matrix of coupling constants in Eq. ( !59|) is Hermitian, e.g. = A|^, then we can 
identify the couplings of Eq. fl2Ul) as 

9s, =gs2 = ^, 9 s, =9s2 = ^ (60) 

for scalar interactions and 

9s, = -9S2 = 9 s, = -9s2 = (61) 

for pseudoscalar interactions. 

To proceed, we need to separate two cases: (i) the lightest FCNC Higgs particle is a 
scalar, and (ii) the lightest FCNC Higgs particle is pseudoscalar. 



1. Light scalar FCNC Higgs 

The case of relatively light scalar Higgs state is quite common, arising most often in Type- 



4i 



42 



3- 



III two-Higgs doublet models (models without natural flavor conservation) 

Bg-Bg Mixing: Given the general formulas of Eq. ( 12T|) . it is easy to compute the contribution 
to AMg^^ of an intermediate scalar (0) with FCNC couplings, 

. ^ W _ 5/|M^J,(C.,m,) (XgV 

^^^^B. - ^ [jf^ j , (6^j 

U(Ci,mb) = ^C3{mb)B3 - (C,{mb)B^ + Crimb)Br) + ^ (C5{mf,)B, + Cs{mb)Bs) , 
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with 'reduced' Wilson coefficients {Ci{fi)} given in Eq. fl22|) . 

—7- fi^fi^ Decay: Comparing Eq. fIBUl) to Eq. fl2^ . we can easily see that the branching 
fraction for the rare decay — t- is zero for the intermediate scalar Higgs, 







(63) 



This is consistent with what was already discussed in Sec. IIII Bl and implies that the FCNC 
Higgs model does not produce a contribution to B^ — )■ yU+yU^ at tree level. The non-zero 
contribution to — )■ fi~^fi~ decay is produced at one-loop level |43|. 



2. Light pseudoscalar FCNC Higgs 



The case of a lightest pseudoscalar Higgs state can occur in the non-minimal supersym- 



metric standard model (NMSSM) 



45 



48| or related models 49|. In NMSSM, a complex 



singlet Higgs is introduced to dynamically solve the /i problem. The resulting pseudoscalar 
can be as light as tens of GeV. This does not mean, however, that it necessarily gives the 
dominant contribution to both B^ — 5^ mixing and the B^ — > fi'^fj^~ decay rate since there 
can be loop contributions from other Higgs states. In the following, we shall work in the 
region of the parameter space where it does. 

B^-B^ Mixing: The contribution to AM^^^ due to intermediate pseudoscalar with flavor- 
changing couplings can be computed using the general formula in Eq. f l2T]) along with the 
identification given in Eq. ( l6T]) . 



48 [mJ ' 

^-C3{m,)Bs + {C,{m,)B, + Cr{m,)Bj) - ^ {C,{m,)B, + Csim,)Bs 



(64) 



fa{Ci,mb) 

with 'reduced' Wilson coefficients Ci{fi) again being defined in Eq. f l22|) . 

B'^ — !■ fi'^^~ Decay: The branching ratio for rare decay can be computed with the help of 
the general formula of Eq. (1241) . 



(a) 



1 PsMl 



327r m^Ep 



1 - 



1/2 



' \D \E ' 
^32 "^22 

M2 



(65) 
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FIG. 2: Branching ratio of Bbo^u+u- as a function of pseudoscalar Higgs mass Ma- Left: Afa = 
1,0.5,0.1 (solid, dashed, dash-dotted lines). Right: X22 = 0.1,0.05,0.01 (solid, dashed, dash-dotted 
lines). In each figure, the yellow shaded area represents excluded parameter space. 

We can now eliminate one of the three unknown parameters (Af^, A^, and Ma) which appear 
in Eqs.( l6^ and (!65|) . We choose to eliminate A^, so 

.(a) _ A ^bA' (i _ ^] (^\ ' reel 

where x^f^ = AM^^^-^/Tb,. As one can see, the unknown factors enter Eq. fIBB]) in the 
combination Afg/^a- It is, however, more convenient to plot the dependence on Ma for 
different values of Afg, which we present in Fig. [2J 

It must be emphasized that the discussion above assumed the absence of large destructive 
interference of the NP and SM contributions to -B^ — mixing. Concrete models where 
such interference is present (and thus the New Physics contribution is larger than the SM 



one) can be constructed 50|. In such models possible contribution to Bs — t- /i^yU could be 
large. 



G. Fourth generation models 

One of the simplest extensions o 
fourth generation of chiral quarks 



the Standard Model involves addition of the sequential 



51 



53l |. denoted for the lack of the better names by t' 



and b'. The addition o: 



the sequential fourth generation of quarks leads to a 4x4 CKM 



quark mixing matrix 5J]. This implies that the parameterization of this matrix requires six 



real parameters and three phases. Besides providing new sources of CP-violation, the two 
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additional 
effects [55 1 



phases can affect the branching ratios considered in this paper due to interference 



There are many existing constraints on the parameters related to the fourth generation 



of quarks. In particular, a fit of precision electroweak data (S and T parameters) 
implies that the masses of the new quarks are strongly constrained to be 59 1 



56 



-l58| 



1 + 



niH 



X 50 GeV, 



V 5 (115 GeV); ^^^^ 
with mf/ > 400 GeV. Here rriH is the SM Higgs mass, which we take for simplicity to be 120 



rrit' — rrib' 



GeV. We also used updated constraints on CKM matrix elements [6] 

The relationship between AMb^ and BB^^fj,+^~ in the model with four generations of 
quarks has been previously studied in detail in [6l|. Here we update their result. The 
branching ratio of Bg — /x^/x^ can be related to the experimentally- measured® xs^ as [61 1 



B 



4m2 \a 



tot\ 
10 I 



(68) 



where the parameter A' is a i^^-mixing loop parameter 6l| . 

A' = r]tSo{xt) + r]t>Rl^So{xt') + 2r]t>Rt>tSo{xt, Xf), (69) 

and Rft = Vt'sV^Ji^/VtsV^l- Bbs can be obtained from Table [11 The definition of the function 
So^Xt^Xt') can be found in Ref. [6l|. The Wilson coefficient C*o* is defined as 

c'^^{^i) = CM + R,tCU^^) (70) 

with Cfg obtained by substituting mf into the SM expression for Cio 62|. The results can 
be found in Fig. [3l As one can see, the resulting branching ratios are still lower than the 
current experimental bound of Eq. ([6l), but for the values of the four-generation CKM matrix 
^bs = of about 0.01, disfavored by 60|, but still favored by [63j, can be quite close 

to it. 



IV. CONCLUSION 

Experiment has determined AM^^ exceedingly well. The Standard Model determination 
provides a consistent value, although with a markedly greater uncertainty (due mainly to 

^ Here we use AA/^^ from Table HI as the separation of NP and SM contributions used in the rest of this 
paper, xb^ — xsm3 + xsma, is not possible due to loops with both t' and c, or u quarks. 
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FIG. 3: Left: branching ratio of Bbo_^,,+, 
values of ' 
(see also 



y^t-^— as a function of the top-prime mass mt> for different 

values of the phase (pt's = 0,11/2, tt (solid, dashed, dash-dotted lines) and A^'^ = ~ 10~^ 6C| 

63]). Right: branching ratio of B^o^s^^+f^- as a function of the CKM parameter combi- 
= and different values of mf = 400 GeV (solid), 500 GeV (dashed), and 
600 GeV (dash-dotted). 



nation A^'^ with 



the dependence on the nonperturbative quantity fsfiss ^"^^ to a lesser extent on the CKM 
mixing element Vts)- We have argued that this fact can be used to constrain NP predictions 
for other processes, such as the Bg — )■ transition considered here. 

We expect this kind of correlation to be a rather general feature of New Physics models, 
provided there is an overlap between the NP parameters which describe AMb^ and (for 
our purposes here) Bg — )■ fi^fi^ . However, given the abundance of New Physics scenarios, 
each with its particular structure, it is not reasonable to expect any universal correlation 
between i?s-niixing and Bg — )■ fi^fi^. Instead, what we have done in this paper is to analyze 
several NP models in detail. In each case, we have first determined the set of unknown NP 
parameters and then, using dynamical assumptions, have been able to reduce (or entirely 
eliminate) the arbitrariness. Analyzing specific NP models this way has two purposes: 
to serve as an instructive example for further study and to see what kinds of numerical 
predictions these particular models yield. 

Not surprisingly, the simplest model (with a single Z' boson) provides a strong correlation 
between AMb^ and Bg — ^ ^^fi^ in which the latter is determined in terms of Mz>- An even 
stronger prediction occurs in the particular version of the Family Symmetry model discussed 
earlier, where a clean determination of Bg — )■ fi'^fi" is obtained. In this instance, a set of 
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reasonable assumptions allows for the initial presence of unknown parameters to be totally 
overcome. A similar, but not quite as fortunate, situation occurs for R-parity violating 
supersymmetry, wherein a reasonable assumption partially reduces the NP parameter set. 
In this case, Bg — )■ fi^fi'^ can be expressed in terms of a ratio of a coupling constant and 
sneutrino mass Mj^. The flavor-changing Higgs model turns out to be less accommodating 
in that no set of assumptions known to us can reduce the original set of three unknown 
parameters. Thus, the constraint from Bg mixing still leaves one with two unknowns (see 
Fig. 12]). We also updated constraints on the models with fourth sequential generation of 
quarks. Of course, additional NP models are available for study, e.g. R-parity conserving 
supersymmetry [25|, and work proceeds on these. 

Finally, as discussed in Sect. Ill, it would be of interest to address the impact of NP 
CP-violating contributions to Bg mixing. Indeed, we plan do so in a future project, but first 
await more accurate data on AF^ or studies of Bg — )■ J /ipcf) transition at LHCb. 
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Appendix A: Choice of the basis and mixing matrix elements 



There are eight Ab = 2 effective operators that can contribute to Sg-mixing. The operator 
basis we shall employ is 



Q3 = ihsR) (bRSL) , Qr = {biSR) {hsR) , 

Qi = (bRSL) (bRSL) , Qs = ihaf^^SR) (bLa^^'sR) 



(Al) 



where quantities enclosed in parentheses are color singlets, e.g. (Pl'-^^sl) = bL,i1fiSL,i- These 
operators are generated at a scale M where the NP is integrated out. A non-trivial operator 
mixing then occurs via renormalization group running of these operators between the heavy 
scale M and the light scale /i at which hadronic matrix elements are computed. 

We need to evaluate the -B°-to-5^ matrix elements of these eight dimension-six basis 
operators. This introduces eight non-perturbative B-parameters {-Bj} that require evaluation 
by means of QCD sum rules or QCD-lattice simulation. We express these in the form 



(Qi) = 




(Q5) 


= flMlB, , 


{Q2) = 




m 


= IflMB, , 


(Q3) = 






= -yiMiB 


{Q4) = 


-yiMB, , 


{Qs) 


= flMlB, , 



(A2) 



7 ; 



where /a. is the Bg meson decay constant and {Qi) = {B^\Qi\B^). 



Ref. 



371] has performed a QCD-lattice determination (quenched approximation) of the 



B-parameters in an operator basis {Oj} which is distinct from the {Qi} of Eq. (lAip . 

01 = %(l + 75)5* 5V(l + 75)5^' , 

02 = F(l + 75)5^ b\l + 75)5^' , O4 = F(l + 75)5^ b\l - 75)5^- , (A3) 

03 = ^(1 + 75)5^' b\l + 75)5^ , O5 = F(l + 75)5^' b\l - 75)5^ . 

Three more operators Oj {i = 6, 7, 8) can be obtained by substituting right-handed chiral 
projection operators with the left-handed ones Oi {i = 1,2,3) in Eq. flA3p . The B^-to-^^ 
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matrix elements of these operators have been parameterized in Ref. 37|] as 



(01) = IMMA ' 

(02) = -IR'JIMIB, , (O4) = 2RlflMlB, , (A4) 

(03) = \RlflMlB, , (O5) = iRlflMlB, . 

Also, the chiral structure of QCD requires that (Oq) = (Oi), (O7) = (O2), and (Og) = (Oa). 

Several of the quantities introduced above are scale dependent, i.e. {-Bj(/i)}, and 
R?s{^i). Throughout this paper, we shall understand all these quantities to be renormahzed 
at a common scale /i = mf, and to simplify notation, we shall denote them simply as {-B^}, 
{Bi} and i?^. In particular, our evaluation at scale fi = nii, of the quantity Rsifi) = 
MBj{mb{n) +ms{n)) yields 

Rl = MlJ (m6(m,) + m,(m,))' = l.ST^Ho , (A5) 



where we have used the input values m{,(m{,) = GeV m and ms{mh) = 0.085 ± 

0.017 GeV 3. 

The two bases {Qi} and {Oj} can be related via Fierz rearrangement, 

01 = 4 Qi , 

04 = 4 Q3 , 

02 = 4 Q4 , (A6) 

05 = -2Q2. 



O3 = -2 ^4 - i Q 



5 ) 



from which we find 





= B, , 


B2 


= iB.R^s 


B3 


— y-D4-flg 


B, 


= B2R'i , 



55 = -ii?f (2^3-55: 



'2 



-Be = -Bi , 

r) 6 rj p2 



(A7) 



Alternatively, the B-parameters can be estimated using the 'modified vacuum saturation' 
(MVS) approach, wherein all matrix elements in Eq. ( 1A2I) are written in terms of (known) 
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List of {B,} 
(in {Qi} Basis) 



Bi = Bq 

B2 

B3 

B4 = B7 
B5 = Bs 



{Bi} from lattice QCD 
(from Ref. [37]) 



0.87 

0.99i?2 
0.80i?2 
0.71/?? 



Bi in MVS 
(from Eq. (jM 



0.87 



0.87 
0.87 



2 p2 



1 I 6 n2 

0.87i?2 
0.87i?? 



TABLE IL Numerical estimates of the B-parameters. The determination from lattice QCD is done 
in MS(NDR). 



matrix elements of {V — A) x (V — A) and {S — P) x (S + P) matrix elements Bb and 



(S) 



(Qi) = ^flMlBs^ 



1 7] 

1 f] 



(Q5) = j;^JImIBb^ rj 
{Qe) = (Qi) , 

{Q7) = m , 



(A8) 



where we take A'^f. = 3 as the number of colors and define 



Ml 



Ri for 5b 



(S) 



B 



B., 



(A9) 



Bb, {rhbirhb) + rhsirhb)) 
It is instructive to compare how well the MVS approximation estimates the recent lattice 
results. We provide such a comparison in Table Ull 
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